The role of the ionic potential in high harmonic generation 
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Recollision processes provide a direct insight into the structure and dynamics of electronic wave- 
functions. However, the strength of the process sets its basic limitations - the interaction couples 
numerous degrees of freedom. In this letter we decouple the basic steps of the process and resolve 
the role of the ionic potential which is at the heart of a broad range of strong field phenomena. 
Specifically we measure high harmonic generation from argon atoms. By manipulating the polar- 
ization of the laser field we resolve the vectorial properties of the interaction. Our study shows that 
the ionic core plays a significant role in all steps of the interaction. In particular, Coulomb focusing 
induces an angular deflection of the electrons before recombination. A complete spatio-spectral 
analysis reveals the influence of the potential on the spatio-temporal properties of the emitted light. 



Electron scattering is one of the main tools to probe 
the structure of matter. It has been successfully used to 
resolve the structure of solids pj, molecular and bi- 
ological 1 3;] complexes. Recent advances in strong field 
light-matter interactions have opened a window into a 
new class of phenomena where electron scattering is ini- 
tiated and manipulated by strong laser fields. In these 
experiments, the strong field leads to tunnel ionization 
which splits the electronic wavefunction into two parts: 
a bound state, associated with the hole left behind the es- 
caping electron, and a continuum state. The free part of 
the electron wavepacket is accelerated in the laser field 
and can recollide with the parent ion. Recollision can 
lead to the emission of optical radiation [3| • This process 
is known as high-order harmonic generation (HHG) [5|. 

In recent years, HHG has been the subject of numer- 
ous studies in which the harmonic spectrum was used to 
retrieve both the spatial [gHlOj and temporal [ll| charac- 
teristics of the electronic states participating in the pro- 
cess. When tunneling occurs from more than one orbital, 
a "hole" wavepacket is launched [12J. Furthermore, a sig- 
nificant population transfer induced by the strong laser 
field may evolve during the optical cycle [ll|. When the 
electron recollides it maps the evolved hole wavefunction 
into the harmonic spectrum. These phenomena are at 
the heart of attosecond science. Their observation, so 
far, has been rather limited albeit they are likely present 
in various complex systems such as large molecules and 
clusters. 

The main challenge is related to the complexity of the 
interaction. All degrees of freedom are coupled by the 
strong field light-matter interaction. How can we isolate 
the required information? One approach is to oversim- 
plify the description. A commonly applied example is 
the Strong Field Approximation (SEA) [14j . which ne- 
glects the influence of the ionic potential not only on 
the free electron trajectory but also on the structure 
of the recolliding wavepacket by treating it as a plane 
wave (plane wave approximation, PWA) . In this letter we 
propose a different approach for probing recollision pro- 



cesses, demonstrating the ability to isolate the role of the 
ionic potential. The ionic potential is in fact involved in 
all stages of HHG. It dictates the structure of the involved 
bound state, defines the tunneling process [15j, and mod- 
ifies both the structure and dynamics of the recolliding 
electron. Isolating its role is a nontrivial challenge since 
the harmonic emission encodes the spatial and temporal 
properties of all the three steps involved in the generating 
process. A direct and unambiguous probing thus relies 
on the following requirements: 1. The study must be per- 
formed with systems in which the bound state remains 
static between ionization and recombination. 2. We have 
to control the recolliding wavefunction separately from 
the bound state. 3. We need to resolve the vectorial 
properties of the induced dipole moment, modified by 
the ionic potential, eliminating the contribution of the 
tunneling probability and the recollision efficiency. 

We address all the above requirements. We perform 
our study with argon atoms, whose ground state can be 
represented by field-free p orbitals even in the presence 
of a strong laser field [l6| . Second, we control the contin- 
uum part of the electronic wavefunction by inducing the 
interaction with an elliptically polarized field. Finally, 
we probe the vectorial properties of the interaction by 
measuring the polarization vector of the harmonics. For 
each harmonic, we separate contributions from two (long 
and short) delayed electron trajectories, which accord- 
ingly, allow us to scan not only the structure but also 
the dynamics of the system. Such separation is achieved 
via a complete spatio-spectral analysis. A semi-classical 
analysis of our measurements reveals the effect of the 
ionic potential on the structure of the recolliding electron 
wavepacket. Furthermore, we show that the longest elec- 
tron trajectories are subject to a dynamical effect from 
the ionic potential, referred to as Coulomb focusing fl7l ]. 

In many high harmonic spectroscopy experiments, the 
angle between the bound state and the recolliding elec- 
tron wavepacket is defined by aligning the considered 
molecules. In order to probe the molecules from differ- 
ent angles, it is possible to change this alignment angle. 
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FIG. 1: (a) High harmonic generation in an elliptical field. 
The red arrow represents the rotating electric field, the green 
fine is the electron trajectory, and the purple arrow indicates 
the polarization direction of the high harmonics, (b) Calcu- 
lated classical short (continuous) and long (dashed) electron 
trajectories leading to emission of harmonic 17 in argon at 
1.2xl0 14 W.cm" 2 and e p = 0.2, in the Strong Field Approx- 
imation. 



However, this is not equivalent to simply changing the 
recollision direction with respect to a fixed bound state: 
rotating the molecule can modify the structure of the tun- 
neling electron wavepacket, the tunneling probability or 
the bound state from which the electrons preferentially 
tunnel (l2| . This problem can be overcome by decou- 
pling the recolliding wavefunction from the bound state 
through the control of its interaction with the laser field. 
Such a scheme has been recently demonstrated using a 
two color orthogonally polarized field Here we use 
an alternative approach, producing high harmonics with 
elliptically polarized light (Fig. [1]). Assuming an adia- 
batic response, tunnel ionization occurs along the instan- 
taneous polarization of the driving laser field. A quanti- 
zation axis aligned in the direction a of the electric field 
at the time of ionization is then selected [HI, [l9[ . The 
elliptically polarized field shifts the freed electron in the 
lateral direction. Due to the initial transverse momentum 
distribution, part of the electron wavepacket can recollide 
with its parent ion. The recollision occurs at an angle (3 
with respect to the main axis of the laser polarization 
ellipse x. Recollision leads to the emission of harmonic 
radiation polarized along the angle x- For spherically 
symmetric bound states the polarization angle is aligned 
along the recollision angle (3 — \- If the bound state is 
not spherically symmetric then the polarization angle x 
is rotated from f3 according to the structure of the con- 
tinuum state and the bound state P3. By controlling 
the ellipticity of the IR field we control the relative angle 
between ionization and recollision. 

The experiments were performed using the 1kHz 800 
nm 35 fs Aurore laser system from CELIA. The ellipticity 
of the generating field is controlled by rotating a zero- 
order half-wave plate in front of a zero-order quarter- 
wave plate. The laser pulses are focused into a 1 kHz 
pulsed argon gas jet. We analyze the produced harmon- 
ics with an XUV flat field spectrometer. We optimize 
the laser focus position with respect to the gas jet to en- 



sure proper spectral and spatial separation of short and 
long trajectories on the detector due to phase matching 
conditions [20[. The polarization state of the high or- 
der harmonics is analyzed by inserting a 45° incidence 
unprotected silver mirror between the grating and the 
detector. The mirror reflectivity is ~ 30 times higher for 
S-polarized radiation than for P. By rotating the rela- 
tive angle between the polarizer and the generating field's 
main axis, we observe a cosine squared modulation of the 
harmonic signal (Malus' law) from which we extract the 
polarization angle of each harmonic [U [22j . 

For each harmonic, the signal from short and long tra- 
jectories is extracted by averaging the polarization angle 
over two separated areas on the detector. The resulting 
polarization angles \ obtained at I « 1.2 x 10 14 W/cm 2 
are shown in Fig. [5] as a function of the laser elliptic- 
ity e p . Several general features can be observed: (i) the 
polarization angle \ generally increases with e p , (ii) the 
polarization angles measured from short and long trajec- 
tories have opposite signs (iii) for the highest harmonic 
orders x remains close to zero for all ellipticities. 

All these features can be qualitatively understood from 
classical analysis. Observation (i) can be easily ex- 
plained: increasing the ellipticity increases the lateral 
displacement leading to larger recollision angles and thus 
to larger \- Feature (ii) can be explained considering the 
classical electron trajectories presented in Fig. Qlb). To 
have a closed trajectory, the electron must be launched 
into the continuum with an initial transverse momentum 
{p y m < 0) counteracting the accumulated momentum 
(p* y leld > 0) . The final transversal momentum is dictated 
by the two components according to: p y ot = p y rn + p^ leld . 
Long trajectories, which interact with almost a complete 
cycle of the laser field, have a vanishing p^ leW , there- 
fore p y ot « p l y m < 0. By contrast, in short trajectories 
pfieid j g s ig n ifi can tly larger than p y m , leading to p y ot > 0. 
Opposite transverse momenta lead to opposite polariza- 
tion angles for l ong and short trajectories as observed 
in (ii) (see also [23[). High harmonics close to the cut- 
off are emitted by coalescing short and long trajectories 
for which p y m is almost totally counterbalanced by p^ leld . 
Accordingly, these trajectories recombine with /3 ~ and 
feature (iii) follows. 

Measuring x provides a direct insight into the vecto- 
rial properties of the interaction. The value of x relies on 
the relative orthogonal recombination dipolc components 
and is intrinsically insensitive to the tunneling probabil- 
ity or the recollision efficiency. Therefore, it allows us to 
isolate the effect of the ionic potential, which we analyze 
hereinafter using a semi-classical theoretical approach. 

Our theoretical model, CTMC-QUEST [24], combines 
a statistical description of electron paths in a Classi- 
cal Trajectory Monte Carlo (CTMC) framework and 
a quantum mechanical description of the recombina- 
tion process using QUantum Electron Scattering Theory 
(QUEST) [25j . In the framework of the single active elec- 
tron approximation, CTMC employs an A-point discrete 
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FIG. 2: High harmonic polarization angles in argon at 
1.2 xlO 14 W.cm" 2 , for short (left) and long (right) tra- 
jectories. Experimental results (black dots) and theoreti- 
cal CTMC-QUEST calculations using scattering continuum 
states (solid blue lines) and plane waves (dot-dashed red 
lines) . 



representation of the electron phase space distribution, 

0(r,p,*) = j? x E^Li <K r - r i(*))<^(P - Pj(t)), in terms 
of N noninteracting trajectories {rj(t), Pj{t)}j=i,..,N- 
Electron dynamics are then monitored by the Liouville 
equation dg/dt — — [g, H] where [ ] is the Poisson bracket 
and H the effective one-electron Hamiltonian which in- 
cludes the model electron-ionic core interaction V c pro- 
posed in [2(| and the laser-atom interaction r ■ E where 

E = (E x ,E y ) 



E (sm(wt), e p coB(wi))/yl + e£. The 

initial conditions are set by defining the distribution 
p(r,p,t — 0) as detailed in [24], to obtain spatial and 
momentum electron densities in close agreement with 
their quantum analog. The classical energy spread of 
the ground state allows mimicking tunneling in terms of 
above barrier transitions [23]. This correspondence has 
been clearly illustrated in stationary barrier cases [28[ 
and quantitatively corroborated in non-stationary ones 
29]. 

Among the N trajectories, we only select those which 
return at time t — t rec to the core with the angle (5. 
Assuming pure selectivity of the bound state achieved by 
the ionization process, as schematically described in Fig. 
1(a), these trajectories contribute to the dipole moment 
according to: 

d u = cos(a) < ^\n\^ Spx > +sin(a) < *kN#3p„ >, 

(1) 

where r» = {x, y}, ^3 Pr are ground-state argon wave- 
functions defined in the laboratory frame, and is 



the continuum wavefunction which describes the return- 
ing electron with momentum k [241 ] . The angle of ion- 
ization, a — tan (Ey(ti on )/ E x (ti on )), determines the 
alignment of the quantization axis with respect to the x- 
axis of the laboratory frame (Fig. 1(a)). Drawing from 
Eq. (1), the HHG polarization angle is finally defined as 
X = Q.btanT 1 {\d x \\d y \ cos<& /[\d x \ 2 - \d y \ 2 }), where $ is 
the phase difference between d x and d y . 

The value of k entering Eq. (1) is directly extracted 
from the returning CTMC trajectory. Evaluating a in 
the classical framework is more critical since a is deter- 
mined by tunnel ionization. Even though it is possible to 
extract approximate ionization times from the classical 



21] 



trajectories |30h33| ]. we rely in practice on the SEA 
to estimate the ij OTi values. 

Figure 3 shows the polarization angles \ calculated 
by CTMC-QUEST, using continuum scattering states 
|24| . The agreement with the measurements is very good. 
Based on this we have taken advantage of our semi- 
classical approach to investigate separately the static and 
dynamical actions of the atomic potential V c on the har- 
monic polarization state. 

First, we investigated the importance of the ionic po- 
tential on the structure of the recolliding electron wave- 
functions. For that purpose we replaced the scattering 
waves 'Fk by plane waves. The results are shown in Fig. 
2 (dashed-dotted lines) . These PWA calculation dramat- 
ically fails in reproducing the observations. 
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FIG. 3: (a) Electron recollision angles at I = 1.2 x 10 14 
W/cm 2 and e p = 0.2, using full (plain lines) and distorted 
(dashed lines) CTMC calculations for short and long trajec- 
tories. The distorted calculations neglect Coulomb focusing, 
(b): Polarization angles for long trajectories of harmonic 15: 
experiment (dots), full (plain line) and distorted (dashed line) 
CTMC-QUEST calculations. 



Second, we studied the dynamical effect of the ionic 
potential on the electron trajectories in the continuum, 
by setting V c = as soon as the trajectories turnaround 
and begin their return to the core. In Fig. 3 we compare 
the resulting recollision angles j3 to those obtained with 
the full CTMC calculation. We observe that short trajec- 
tories are quite insensitive to the potential. By contrast, 
V c induces a significant decrease of the recombination an- 
gles associated with long trajectories. As we go to lower 
harmonic orders, the decrease in /3 becomes more signif- 
icant (up to 40% for the 13 th harmonic). This angular 
shift is a signature of Coulomb focusing [17], an effect 
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FIG. 4: Spatio-spectral distributions of intensity and po- 
larization angle in the far-field, for harmonic 15 driven by 
e p — 0.22 . Short trajectories are collimated and long trajec- 
tories appear as a ring. 



which has been suggested as an explaination of enhanced 
double ionization. To our knowledge, our results are the 
first demonstration of the effect of core focusing on HHG. 

The different importance of Coulomb focusing to short 
and long trajectories can be intuitively explained. Short 
trajectories return to the ionic core under the combined 
influence of the laser and ionic fields, which both accel- 
erate them towards r = 0. In contrast, long trajectories 
are slowed down since the two fields counteract each other 
before recombination [32l - i3~4T | . Since cutting off a coun- 
teracting action has more consequences than canceling a 
concurrent one, V c has more influence on the returning 
dynamics of the long trajectories than of the short ones. 

Whereas Coulomb focusing has no effect on the po- 
larization angle for short trajectories, it clearly reduces 
the polarization angle associated with the longest elec- 
tron trajectories (Fig. |3]Jb)). From a constructive point 
of view, this shows that the SFA can be safely used 
to describe above-threshold HHG associated with short 
trajectories, provided that accurate scattering states are 
used in the recombination step to overcome the above- 
mentioned PWA shortcomings. On the other hand, care 
must be taken when employing SFA to long trajectories 
even though accurate continuum states are used. 

As in most high harmonic spectroscopy experiments, 
we have so far approximated the harmonic signal result- 
ing from spatial and temporal averaging by the single 
atom response at a given effective laser intensity. The 
excellent agreement between our simulations and the ex- 
periment confirms the validity of this approach. How- 
ever, we can take one step further by analyzing the en- 
tire spatio-spectral distribution of the polarization an- 
gle in our experiment. Figure [4] shows this distribution 
for harmonic 15 as it appears on our detector. A clear 
separation between short (central part) and long (ring) 



trajectories is observed [20[. Measuring the polarization 
state of each point on our detector, we resolve the com- 
plete spatio-spectral distribution. This analysis confirms 
the dramatic difference in polarization states of short and 
long trajectories. In addition, it provides a direct visu- 
alization of the intensity-dependence of the polarization 
angle. In our experiment the field intensity changes both 
temporally and spatially. As a result, the generation of 
each harmonic order is accompanied by a spatio-temporal 
chirp. The temporal chirp provides a simple mapping 
between time and frequency: the spectral maximum cor- 
responds to the maximal temporal intensity, whereas the 
spectral wings correspond to the leading and falling edges 
of the pulse. A similar mapping occurs between the near 
field and far field's spatial profiles. Our measurements 
show that for each trajectory, the polarization angle de- 
creases as we move away from the center of the spatio- 
spectral distribution, i.e. as the laser intensity decreases. 
This is consistent with the evolution of the recombination 
angle as we go from the plateau to the cutoff region (Fig. 
[2]). Such an analysis can be easily applied in a broad 
range of experiments, resolving the spatio-spectral re- 
sponse of the interaction to a variable parameter - pump- 
probe delay, two-color phase or carrier-envelop phase. 

The simultaneous spectral and spatial resolution of the 
detection allows us to measure both short and long tra- 
jectories, which can probe a system at a given electron 
energy from two different recollision angles, associated 
with two different recollision times. Such a measurement 
naturally enables the decoupling between structural and 
laser-induced information in high harmonic spectra. In 
this respect, our study contributes to HHG-based spec- 
troscopic methods that aim at resolving the electronic 
structure and dynamics of atomic and molecular systems. 

Although we have focused on argon to reveal the role of 
a static ionic core, this work can straightforwardly be ex- 
tended to more complex systems. One important aspect 
will then be the sensitivity to ionization times through 
the direction of the quantization axis a, which should 
enable studying the influence of core polarization and 
Stark shifts on tunneling [l5j]. This scheme can also be 
extended to resolve the influence of the ionic potential in 
molecular systems. Rotating the molecular axis will pro- 
vide an additional information about the 2D properties 
of the potential. 
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